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Abstract 

We present in this paper a rigorous method to derive the nonlinear 
Fokker-Planck (FP) equation of anomalous diffusion directly from a 
generalization of the principle of least action of Maupertuis proposed 
by Wang [1] for smooth or quasi-smooth irregular dynamics evolving in 
Markovian process. The FP equation obtained may take two different 
but equivalent forms. It was also found that the diffusion constant 
may depend on both q (the index of Tsallis entropy [2]) and the time 
t. 

PACS : 02.50.-r; 05.20.-y; 05.70.-a 

The Fokker-Planck equation is a differential equation describing the time 
evolution of probability distribution of state during stochastic processes. 
The FP equation and its generalizations play very crucial roles in statis- 
tical physics. The FP equation is not only applicable to the systems near the 
thermal equilibrium, but to the systems far from the thermal equilibrium as 
well. This latter application has special meaning in dealing with a large class 
of self-organized, complex dynamical systems. In this sense, the FP equation 
not only describes stationary properties but also the dynamics of evolving 
systems. 
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The FP equation was first derived by Fokker [3] and [4] as one to describe 
Brownian motion. Later on, many books and review articles were published 
[5, 6, 7, 8, 9]. The usual way of deriving the FP equation starts from calculat- 
ing the transition probability P(x,t + r\x',t) for small r, where the particle 
travels from x to x' . There are various ways to derive the expression of the 
transition probability [10]. In this paper, we employ the concept of max- 
imum path information, related to non-extensive Tsallis entropy, to derive 
the expression of the transition probability for the motion of particle under 
the influence of external forces. On the basis of the transition probability, 
an nonlinear FP equation can be obtained. 

A path information based on Shannon entropy [11] has been defined as 

[1] 

w 

H s (a,b) = -^2p ab (i)\np ab (i), (1) 
i=i 

where Pab(i) is called the transition probability that a system moving from 
point a to point b will choose path i among all possibly existed w paths. 
In the case of non-extensive statistics, the corresponding path information 
naturally adopts the formula of Tsallis entropy [2] 

H t (a, b ) = -k± Mt) - p ' J: '\ (2) 

1=1 1 1 

where q, the entropy index, specifies a particular statistics. In general, the 
larger the path information, the less we know about paths states of the 
system. 

Of course, the transition probability p a b{i) satisfies the following normal- 
ization condition 

to 

5>*(0 = i- ( 3 ) 

1=1 

For classical dynamical systems we also suppose each possible path is char- 
acterized by its action A ab {i) 

A ab (i) = f Li(t)dt, (4) 

J tab 

where Li(t) is the Lagrangian of the system at time t via the path i. The 
average action is represented by 

w 

(A ab ) = J2 A ab(i)Pab(i). (5) 

1=1 
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In order to obtain the form of path probability, we seek to optimize the 
path information H t (a, b) under the constraints of Eqs. (3) and (5). That is, 

w w 

5(H t (a,b) + aY,Pab(i) + vJ2 A *b(i)Pab(i)) = (6) 

i=l i=l 

Through a simple algebra, the optimization yields the following expression 
of path probability 

Pab (i) = ^[l-(l-q)vA ab ( i )]^, (7) 

where Z q = £^[1 - (1 ~ q)r]A ab {i)]^ . 

In order to obtain a general derivation of FP equation at the existence 
of any form of external forces (drifts), we adopt here the Euler's method to 
calculate the action. The detailed method is as follows. The path through 
which the particle travels from point a to point b is cut into TV segments with 
each having a spatial length Axk = (xk — Xk-i) (k = 1, 2, N). t = tf.-tf.-i 
is the time interval spent by the system on every segment. According to the 
theorem of large numbers, the fluctuation of calculation will go to as N 
approaches infinity. The action on the segment k is simply 

m(Ax fc ) 2 Ax k 

2t ~2T k ~ U \ Xk - 1 ) 1 ^ ( 8 ) 

where F k = — (f^)fc and U(xk-i) is the potential energy at the point x k -\- 
Here in this paper F k and U(xk-i) will be considered as constant. From now 
on, we will write U(xk-i) as U for simplicity. 

By using Eq. (7) the transition probability Pk/k-i from point k — 1 to 
point k via the path % can be written as 

Pk/k ~ lil) = z g (M-i) {1 - (1 -^ [ ^T- 

+^Ax k -Ut]}^, (9) 

where Z q (k, k—1) can be calculated from the normalization condition p k / k -i(i)dxk 
1 

r m(Aa; fc ) 2 

-oo 2t 



/+oo j 
dx k {l - (1 - q)v[- 
-oo 



+^A Xk -Ut}}^. (10) 
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Introducing the methods by Tsallis and Prato [12, 13], after a tedious 
calculation, we obtain the exact form of Z q (k, k — 1) 



Z q (k, k-l) = A(q)^—[1 -(q- l) V {-^- + Ut)]£*, (11) 



where A(q) can be written as 



A(q) 



It is not difficult to prove that Z q (k,k — 1) restores to Zi(fc, k — 1) at the 
q — > 1 limit, which is 



Zi(fc, fc - 1) = e*pfo(-^- + Ut)]J — . (12) 

8m y rjm 

Hence the transition probability pk/k-iif) has the form 

Pfc/fc _i(i) = B(g)r 1/2 [i-(?-iM|^ + ^#{i- 

8m 

(l- ? ),[^^ + fA, t -«]}A, (13) 

where 5(g) = A'^q)^. 

Now we are ready to derive the FP equation for the system travelling 
through the fc-th segment of path i connecting points a and b. It is readily 
that 

<% ^ 7 U A 2 1 2 v 8m 

tO + A^ v [^(Ax k y - ^ + U}}, (14) 

where A x = l-( q -l) v (^ + Ut), and A 2 = l-(l- g ) r/ [|(Ax fc ) 2 + ^if^t- 
Ut]. We also have 

= -p^m^mr'Ax, + ^), (15) 
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and 
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x^M^Ax. + ^f], (16) 

where 7 is a constant that might depend on q. 

Combining the equations (14) and (15), one obtains the following expres- 
sion 

+ F k -^—)p k / k _i{i) = ' [ui(x k ,t,q) + v(x k ,t,q)], (17) 

where 

Fht 

Mi(x fc ,t,g) = mt~ l -rjA^imt^Axk + -^-f (18) 



and 



3F 2 t 2 FH 
v(x k ,t,q) = (—^ + U)(3-q)r ] mAi l + 2 V mA2 l (^- 
8m 2 

+ Sl + ^F ! A Eh + FtAxk _ u) ( 
8m t 

Writing (16) in another form one gets 

d2[Pk Q x 2 {l)V = -7^W-i(*)MM,<?,7), (20) 



where 



Fit 

u 2 (x k , t, q, 7) = mt~ l - (7 - 1 + q) V A^ 1 (mr 1 Ax k + — f. (21) 



It is obvious that u 2 (x k , t, q, 2 — q) = Ui(x k , t, q). 

Relating Eqs.(17) and (20), together with u 2 (x k , t, q, 2 — q) = Ui(x k ,t,q), 
one obtains the following equation 



{ di + Fk dx^ )Pk/k - l{t) = D{q > t) dxl ' (22) 

where 



D fat) = J—B^Af- q)/2 6 l -^ 2 [l - -^hM ]. (23 ) 
1 ' 2 V m 1 1 u 1 (x k ,t,q)(2-q) i 1 ; 



One can check that D(l,t) = l/2r]m, which is consistent with the results 
in [1]. Apparently, Eq. (22) is the exact FP equation for the system in an 
infinitesimally interval in the existence of external forces. 

Besides Eq. (22), the FP equation can also take another form, 



d „_„ d 

x 
of 2 



^IPk/k-iii)?" 1 = -^{F k \pk/k-i{i)] 2 - q )} + D'{q,t) 



xj^lPk/k-iii)] 2 " 1 , (24) 



where 



= ■ ]■ (25) 

2r]m u 1 (x k ,t,q)(2 - q) 

We note that D(q,t) in Eq. (22) and D'(q,t) in Eq. (24) are both q 
and t dependent. The dependence on q is a direct consequence of the non- 
extensive statistics where q is the identity of the system described. It has 
been shown above that when q = 1, the normal diffusion constant can be 
restored. The dependence on t is also quite natural because we are now 
dealing with evolutionary processes where the phase space through which 
the diffusion occurs is changing with time. As f -> oo, one readily obtains 
the diffusion constant for the stationary state. 

The nonlinear FP equation derived above, Eq. (22) and Eq. (24), is well 
applied to describing the evolutionary processes and stochastic processes of a 
large class of self-organized systems that are far from thermal equilibrium, as 
well as chemical equilibrium, such as transportation and diffusion occurred 
in fractal or curved space. For example, it can be employed to describe 
the broad range of markets and exchanges characterized by the anomalous 
(super) diffusion and power- law distributions [15]. Another hope is that this 
equation can be applied to the complex biological systems where evolution 
and anomalous diffusion are taking place from time to time. Compared to 
the normal FP equation and some of its other nonlinear forms [16], our FP 
equation is more general because it can describe both regular dynamics and 
irregular dynamics that occurred in a large category of non-equilibrium and 
chaotic systems [17, 18]. Another important feature of our FP equation is 
that the diffusion coefficient is both q and t dependent. 
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